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FOREWORD 


The  finite  element  analysis  of  the  pile  foundations  of  the  Colum- 
bia Lock  was  sponsored  by  the  U.  S.  Army  Engineer  District,  Vicksburg 
(VED),  and  was  conducted  by  the  Soils  and  Pavements  Laboratory  (S&PL), 

U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES),  Vicksburg, 
Mississippi. 

The  investigations  described  herein  were  performed  by  Dr.  C.  S. 
Desai  and  Dr.  L.  D.  Johnson,  Research  Engineers,  Soil  Mechanics  Division 
(SMD).  Mr.  C.  M.  Hargett,  Engineering  Division,  VED,  monitored  the 
project  and  provided  valuable  comments  and  advice.  Mr.  W.  C. 

Sherman,  Jr.,  Research  Consultant,  SMD,  and  Mr.  S.  J.  Johnson,  Special 
Assistant  to  the  Chief,  S&PL,  provided  a technical  review  of  the  report. 
The  report  was  reviewed  and  approved  by  the  Vicksburg  District  prior  to 
publication.  Mr.  V.  M.  Agostinelli,  VED,  provided  useful  comments. 
During  the  investigation,  Mr.  J.  P.  Sale  was  Chief,  S&PL;  Mr.  C.  L. 
McAnear  was  Chief,  SMD. 

Directors  of  WES  during  the  analysis  and  preparation  and  publica- 
tion of  this  report  were  BG  E.  D.  Peixotto,  CE,  and  COL  G.  H.  Hilt,  CE. 
Mr.  F.  R.  Brown  was  Technical  Director. 
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NOTATION 

A Projected  area  of  pile 

A Equivalent  area  of  piles  battered  inward 

Ago  Equivalent  area  of  piles  battered  outward 

c Mohr-Coulomb  cohesive  strength  parameter 

c Adhesion 

a 

d Experimentally  determined  parameter 

D Depth  of  overburden 

E Modulus  of  elasticity 

E Equivalent  modulus  of  elasticity  for  piles  battered  inward 

Ego  Equivalent  modulus  of  elasticity  for  piles  battered  outward 

E.  Initial  modulus 

1 

E Tangent  elastic  modulus 

F Experimentally  determined  parameter 

G Experimentally  determined  parameter 

k^  Initial  tangent  stiffness 

k Normal  stiffness 

n 

k Residual  stiffness 

r 

k^  Tangent  shear  stiffness 

K Hyperbolic  loading  parameter 
Kj  Experimentally  determined  parameter  for  interface 

Coefficient  of  earth  pressure  in  compression 
Ku  Hyperbolic  unloading  parameter 
L Length  of  piles 

L ^ Equivalent  length  of  piles  battered  inward 

Equivalent  length  of  piles  battered  outward 

m Number  of  rows  in  a monolith 
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Experimentally  determined  parameter;  also,  total  number  of 
piles 

Number  of  piles  battered  inward 
Number  of  piles  battered  outward 
Atmospheric  pressure 
Failure  ratio 

Equivalent  stiffness  per  strip  for  piles  battered  inward 

Equivalent  stiffness  per  strip  for  piles  battered  outward 

Stiffness  per  strip  for  piles  battered  inward 

Stiffness  per  strip  for  piles  battered  outward 

Total  axial  stiffness 

Unit  weight 

Unit  weight  of  water 

Angle  of  friction  between  lock  or  pile  material  and  surround- 
ing soils 

Parameter 

Parameter 

Poisson's  ratio 

Poisson's  ratio  at  failure 

Tangent  Poisson's  ratio 

Major  principal  stress 

Confining  pressure 

Normal  stress 

Horizontal  stress 

Vertical  stress 

Stress  difference 

Shear  stress 

Mohr-Coulomb  (angle  of  internal  friction)  parameter 


CONVERSION  FACTORS,  BRITISH  TO  METRIC  UNITS  OF  MEASUREMENT 

British  units  of  measurement  used  in  this  report  can  he  converted  to 
metric  units  as  follows : 


Multiply 


feet 

miles  (U.  S.  statute) 
square  feet 

pounds  per  square  foot 
pounds  per  cubic  foot 
tons 
kips 


JiL 


0.301*8 
1.6093UU 
0.09  290  301* 
1*7.88026 
16. 0181*6 
907.181*7 
1*1*1+8.222 


To  Obtain 


meters 
kilometers 
square  meters 
pascals 

kilograms  per  cubic  meter 

kilograms 

newtons 
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SUMMARY 

The  Columbia  Lock  was  designed  as  a gravity-type  structure  in 
which  the  load  is  transferred  essentially  through  the  foundation  piles. 
Field  performance  of  the  lock  pile  system  was  measured  in  terms  of 
settlements  and  loads  in  the  piles.  The  structure  was  analyzed  using 
Hrennikoff's  method,  which  is  based  on  a number  of  simplifying  assump- 
tions such  as  linear  material  behavior  and  linear  distribution  of  loads 
in  the  pile  groups.  The  observations  in  terms  of  distribution  of  loads 
in  the  piles  did  not  show  good  correlation  with  the  distribution  pre- 
dicted by  Hrennikoff's  method. 

In  order  to  evolve  a more  rational  procedure  for  design  analysis 
and  for  improved  understanding  of  behavior  of  the  pile  foundations,  it 
was  proposed  to  analyze  the  system  using  the  finite  element  method.  The 
method  simulated  major  steps  of  the  construction  sequence  including  in 
situ  stresses,  dewatering,  excavation,  construction  of  piles  and  lock, 
backfilling,  filling  the  lock  with  water,  and  development  of  uplift 
pressures.  Nonlinear  behavior  of  soils  and  of  interfaces  between  the 
lock  and  surrounding  soils  and  between  piles  and  the  foundation  soils 
was  introduced  into  the  analysis. 

For  simplicity  and  to  avoid  undue  amounts  of  manpower  and  computer 
efforts,  the  three-dimensional  system  was  idealized  as  a structurally 
equivalent  two-dimensional  plane  strain  system. 

The  numerical  results  in  terms  of  loads  in  the  piles  and  settle- 
ments showed  satisfactory  correlation  with  the  observations.  It  is 
believed  that  the  approximate  two-dimensional  analysis  can  provide  solu- 
tions of  practical  accuracy. 

An  available  computer  code  which  was  developed  for  U- frame  locks 
was  modified  to  handle  the  gravity-type  structures,  particularly  the 
foundation  piles  and  interfaces.  Although  the  results  were  found  to  be 
satisfactory,  it  seems  that  the  code  will  need  further  modifications. 
This  will  involve  proper  introduction  of  interfaces  and  their  properties 
and  a scheme  that  can  avoid  excessive  deformations  of  "air"  elements 
that  need  to  be  introduced  during  various  excavation  and  construction 
sequences . It  appears  that  the  existence  of  air  elements  renders  the 
total  stiffness  matrix  ill-conditioned  and  hence  the  excessive  deforma- 
tions. It  may  be  possible  to  introduce  a scheme  such  that  the  need  for 
air  elements  is  avoided  in  most  computations.  Moreover,  the  work  can  be 
continued  further  to  analyze  the  interaction  effects  and  the  behavior  of 
the  structure  itself. 
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FINITE  ELEMENT  ANALYSIS  OF  THE  COLUMBIA 


LOCK  PILE  FOUNDATION  SYSTEM 


PART  I:  INTRODUCTION 


1.  Many  lock  walls  are  designed  as  gravity-type  structures  sup- 
ported on  battered  pile  groups.  A system  incorporating  piles  with  al- 
ternate rows  battered  in  opposite  directions  poses  a three-dimensional 
problem.  Although  three-dimensional  formulations  and  codes  based  on 
the  finite  element  (FE)  method  are  possible,  their  use  for  the  lock 
pile  foundation  problem  would  involve  formidable  amounts  of  human  and 
computer  efforts.  The  three-dimensional  problem  is  therefore  replaced 
by  a structurally  equivalent  two-dimensional  system.  To  the  authors' 
knowledge,  a problem  of  this  type  has  not  been  solved  previously  using  a 
mamerical  (FE)  method. 

2.  The  Columbia  Lock  walls  were  designed  as  a gravity- type  struc- 

12  3 

ture.  ’ Hrennikoff's  method  was  used  for  design  analysis  and  to  eval- 
uate distribution  of  loads  in  the  pile  groups . Field  observations  dur- 
ing and  subsequent  to  the  construction  of  the  lock  were  not  in  agreement 
with  the  distribution  of  loads  in  the  pile  groups  predicted  by  Hrenni- 
koff's method.  It  was,  therefore,  proposed  to  use  the  FE  method  to 
simulate  as  closely  as  possible  various  sequences  of  construction  and  to 
predict  histories  of  settlements  and  distribution  of  loads  in  the  pile 
groups.  The  sequences  simulated  in  the  analysis  included  dewatering, 
excavation,  installation  of  piles,  construction  of  the  piles  and  lock, 
backfilling,  filling  the  lock  with  water,  and  development  of  uplift 
pressures.  These  sequences  were  preceded  by  computation  of  in  situ 
stresses  under  gravity  loading.  The  numerical  predictions  are  compared 
with  observations,  and  conclusions  regarding  design  analysis  are 
presented. 
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PART  II:  DETAILS  OF  LOCK  AND  FOUNDATION  SOILS 


3.  Columbia  Lock  is  located  in  a cutoff  between  miles  131.5*  and 

13L.5  on  the  Ouachita  River  near  Columbia,  Louisiana.  Fig.  1 shows  a 

typical  section  through  the  lock  chamber  togetner  with  foundation  and 
2 1+  5 

surrounding  soils.  ’ Details  of  two  monoliths,  10-L  and  11-R,  that 
were  instrumented  and  are  considered  in  the  FE  analysis  are  shown  in 
fig.  2.  This  figure  also  shows  details  of  the  strain  gages  installed 
on  the  steel  H-piles.  A plan  of  the  lock  and  details  of  the  instru- 
mented monoliths  are  shown  in  fig.  3. 

Subsoils 


L.  The  subsoils  beneath  the  lock  consist  essentially  of  cohesive 

backswamp  deposits  and/or  cohesionless  substratum  deposits  beneath  the 

east  wall  and  Tertiary  deposits  interfingerea  with  colluvium  and  sub- 

U 6 7 

stratum  deposits  beneath  the  west  wall,  fig.  1.  * ’ The  backswamp  de- 
posits consist  of  stiff  to  very  stiff  fat  clays  and  very  soft  to  medium 
silty  clays  and  silty  sands.  Substratum  deposits  consist  of  dense  to 
very  dense  sands.  The  Tertiary  soils  can  be  divided  into  two  parts. 

The  upper  part,  to  el  -20.00,**  consists  of  dense,  gray,  massive  to 
poorly  bedded  silts,  silty  sands,  and  fine  sands  with  scattered  horizon- 
tal layers  of  stiff,  gray,  and  grayish-brown  clay.  The  lower  portion  is 
composed  of  thicker  zones  of  stiff,  gray  to  brown  clay  mixed  with 
thinner  layers  of  dense  to  very  dense,  gray,  fine,  and  silty  sands. 


Stress-Strain  Parameters 


Soil  behavior 

5.  Results  from  a number  of  triaxial  tests  performed  previously 


* A table  of  factors  for  converting  British  units  of  measurement  to 
metric  units  is  presented  on  page  ix. 

**  All  elevations  (el)  cited  herein  are  in  feet  referred  to  mean  sea 
level. 
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Details  of  instrumented  monoliths 


1 2 

on  the  Tertiary  clays  were  available.  ’ For  the  FE  analysis,  the 

stress-strain  curves  were  simulated  by  using  the  hyperbolic  formulation. 

Details  of  the  hyperbolic  formulation  are  given  in  various  publica- 
ft— 1 s 

tions°”  5 and  are  not  repeated  herein.  According  to  this  concept,  the 
tangent  elastic  modulus,  E , and  the  tangent  Poisson's  ratio,  v , 

L»  X. 

are  given  by 


Et‘(1-xi)2Ei  (1) 


in  which 

E.  = initial  modulus 

l 

K = hyperbolic  loading  parameter 
Pa  = atmospheric  pressure 
= confining  pressure 
n = exponential  parameter 
= stress  difference 

G,  F,  and  d = parameters  determined  from  laboratory  tests 

and 

R (a  - a )(l  - sin  0) 

i = _i — ± 2 (■? 

1 2c  cos  0 + 2a ^ sin  0 

Here  R^,  = failure  ratio  and  c and  0 = Mohr-Coulomb  strength 
parameters . 

6.  The  hyperbolic  parameters  for  other  soils  were  adopted  from 
test  data  on  similar  soils  reported  in  previous  studies. Table  1 
shows  the  hyperbolic  parameters  and  other  properties  of  soils  adopted 
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for  the  FE  analysis.  The  values  of  Kq  , the  coefficient  of  earth  pres- 
sure in  compression,  for  clay  soils  were  estimated  on  the  basis  of 

values  for  similar  overconsolidated  soils.  The  question  of  field  evalu- 
c 

ations  of  Kq  is  wider  in  scope  and  will  require  considerable  research 
in  the  future.  Since  adequate  information  on  volumetric  behavior  was 
not  available,  v , the  Poisson's  ratio,  was  assumed  to  remain  constant. 
Moreover,  the  analysis  does  not  permit  values  of  v >_  0.5  , and  it  is 
believed  that  the  gross  predicted  behavior  may  not  be  affected  signifi- 
cantly if  v were  treated  as  variable.  Values  of  the  Young's  modulus 

8 2 

E for  concrete  were  computed  as  W 10  lb/ft  , and  the  equivalent  E 

7 

for  piles  in  monoliths  11-R  and  10-L  was  computed  as  6.8  * 10  and 

11.8  x io^  lb/ft2  for  piles  battered  out,  and  7.5  x 10^  and  13.3  x 10^ 

2 

lb/ft  for  piles  battered  in,  respectively.  The  densities  of  equivalent 
piles  (6  and  7 in  table  l)  were  computed  by  equating  the  weights  of  the 
equivalent  piles  and  the  H-piles.  Computation  of  the  equivalent  E for 
piles  is  given  subsequently.  The  value  of  E for  air  elements  was  set 
at  10  lb/ft2. 

Interface  behavior 

8 9 11  13 

7.  Interface  elements0’  ’ ’ were  introduced  between  the  lock 

and  surrounding  soils  and  between  the  equivalent  piles  and  foundation 
soils.  The  FE  mesh  for  the  lock  and  foundation  is  shown  in  fig.  U; 
part  of  the  mesh  with  details  of  interface  elements  is  shown  in  fig.  5. 
Details  of  the  soil  layers  in  fig.  U were  obtained  from  boring  logs  in 
the  neighborhood  of  monoliths  10-L  and  11-R. 

8.  The  nonlinear  interface  behavior  was  simulated  by  using  the 
hyperbolic  parameters  adopted  from  direct  shear  test  results  for  simi- 
lar soils  (sands  and  clays)  with  pile  material  (concrete  and  steel)  at 


other  sites. 
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kst  ■ (*  - h)\ 


(h) 


where 


/ v r 
/o  \ 

i , = K.y  (—  ) 

1 J W \Pa  / 


_ V 

and  ^2  c + a tan  6 
a n 
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Details  of  piles  and  interfaces 


and 


kst  = tangent  shear  stiffness 

k.  = initial  stiffness 
1 

Y = unit  weight  of  water 
w 

a = normal  stress 
n 

t = shear  stress 

c = adhesion 
a 

6 = angle  of  friction  between  lock  or  pile  material  and  sur- 
rounding soils 


Kj  , n , and  are  determined  from  test  data.  Hyperbolic  parameters 

for  the  interfaces  are  shown  in  the  following  tabulation.  As  soon  as 

the  shear  stress  at  the  interfaces  equalled  or  exceeded  the  Mohr-Coulomb 

strength,  the  value  of  shear  stiffness  was  set  equal  to  the  residual 

stiffness,  k (see  tabulation).  The  normal  stiffness,  k , was 
r 8 3 n 

assigned  a high  value  of  10  lb/ft-3  before  failure  and  a small  value 

3 

of  10  lb/ft  after  failure.  Interfaces  between  the  same  material  were 
kept  inoperational  by  assigning  both  k and  k^  very  high  values  of 

S X 


108  lb/ ft3. 


Material 

kr 

lb/ft3 

6_ 

c 

a 

K. 

J 

n 

fif 

Remarks 

12 

10 

33 

0 

7.5  x 

io4 

1.0 

0.87 

Lock  backfill 

13 

10 

26 

0 

2.5  x 

io4 

1.0 

0.87 

Piles,  sand 
substratum 

lit 

10 

20 

0 

ro 

VJl 

X 

io4 

1.0 

0.87 

Piles,  Tertiary 
clay 

10 


where 


A = projected  area  of  pile 
n = total  number  of  piles 

= number  of  piles  battered  inward 
nQ  = number  of  piles  battered  outward 
A = projected  area  of  pile 

E = modulus  of  elasticity  for  the  steel  H-piles 
L = length  of  pile  (all  piles  are  of  equal  length  and  area) 

For  the  analysis  herein,  the  length  of  the  piles  was  assumed  to  be 
70.0  ft.  It  is  noted,  however,  that  the  magnitudes  of  the  length  can 
be  different  if  the  piles  are  assumed  to  be  point-bearing  or  frictional. 
If  the  total  stiffness  is  assigned  to  m number  of  strips  or  rows  in  a 
monolith,  fig.  6,  the  stiffness  per  strip  for  piles  battered  inward  and 
outward  is 


s . 
1 


s 

o 


^i_  AE 
m L 

% AE 
m L 


(6a) 

(6b) 


where 

s^  = stiffness  per  strip  for  piles  battered  inward 
sq  = stiffness  per  strip  for  piles  battered  outward 
The  equivalent  stiffness  of  a strip  is  defined  as 

A .E  . 
ei  ei 

Sei  L . 

ei 


s 

eo 


A 


eo  eo 
L 


eo 


(7a) 


(7b) 


where 


s . 
ei 

A . 
ei 

E . 
ei 

L . 
ei 


equivalent 

equivalent 

equivalent 

equivalent 


stiffness  per  strip  for  piles  battered  inward 
area  of  piles  battered  inward 

modulus  of  elasticity  for  piles  battered  inward 
length  of  piles  battered  inward 
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s = equivalent  stiffness  per  strip  for  piles  battered  outward 
A^o  = equivalent  area  of  piles  battered  outward 

= equivalent  modulus  of  elasticity  for  piles  battered  outward 
Lgo  = equivalent  length  of  piles  battered  outward 

11.  Two  FE  analyses  were  performed,  one  with  piles  battered  out- 
ward and  one  with  piles  battered  inward;  and  properties  of  each  were 
evaluated.  In  both  cases  a unit  width  of  strip  was  chosen  for  conve- 
nience, and  equivalent  properties  for  unit  length  were  evaluated.  The 
following  are  the  computations  for  the  two  cases,  fig.  7a  and  b. 

12.  The  cross-sectional  area  of  the  H-pile  used  in  the  foundation 

1 2 

was  0.ll*9  sq  ft.  ’ In  the  case  of  piles  battered  inward,  there  were 
six  rows  of  piles  in  monolith  10-L  and  five  rows  in  monolith  11-R.  For 
monolith  10-L  piles  with  inward  batter,  shown  in  fig.  7a,  n^  = 53  , 

A = 0.1U9  sq  ft,  E = L.Ol*  x 109  lb/ft2,  m = 6 , A = 1 x 1+0  (unit 
width  multiplied  by  length  of  the  monolith)  = 1*0  sq  ft,  and  L = L . 


Therefore,  equating  equations  6a  and  7a, 


ei 


ei 


53  x 0.11*9  x k.Ok  x io9 

6xl+o  ‘ 1* 


33  x io8  lb/ft2 


(8a) 


Similarly,  for  monolith  11-R  piles  with  inward  batter 


25  x 0.1U9  x l+.OU  x 10”  7 , 2 

Ee.  = ^ x-.^u HL.  = 7.5  x 10  f lb/ft* 


(8b) 


For  piles  with  outward  batter,  fig.  7b,  the  equivalent  moduli  for  mono- 
liths 10-L  and  11-R  are,  respectively: 


and 


. 38  » O.U?  x hM  « 10*  . la8  , 108  2 

eo  5 x 1*0 


= 2Li^Ui!2^2!LJLWl , 6.6  x 107  lb/ft; 
eo  6xl+o 


(9a) 


(9b) 
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PART  IV:  FINITE  ELEMENT  ANALYSIS  AND  RESULTS 


13.  Finite  element  analyses  of  U-frame  locks  were  performed  by 

8 8 
Clough  and  Duncan.  The  computer  code  developed  previously  was  uti- 
lized for  the  present  analysis.  However,  a number  of  modifications  and 
corrections  were  introduced.  The  stresses  in  the  interface  elements 
were  made  to  be  consistent  with  those  in  the  surrounding  soil  elements. 
Modifications  were  made  to  introduce  interfaces  and  piles  in  the  foun- 
dation such  that  proper  interfaces  became  operational  as  soon  as  the 
buildup  sequence  was  initiated.  Similar  changes  were  made  for  inter- 
faces that  became  operational  as  the  structural  buildup  and  backfill 
proceeded. 

Simulation  of  Construction  Sequences 

ll+.  The  steps  or  cases  in  the  simulation  of  construction  se- 
quences are  shown  in  table  2.  A cycle  of  the  finite  element  procedure 
was  first  performed  to  compute  initial  in  situ  stresses  caused  by  the 
gravity  loading.  Dewatering  was  then  simulated  to  lower  the  water 
table,  which  was  followed  by  excavation  in  three  steps.  The  materials 
of  elements  occupying  the  locations  of  battered  piles  were  subsequently 
changed  from  soil  to  steel  piles  with  equivalent  properties.  The  lock 
structure  was  then  constructed  in  three  steps.  Sequences  11  and  12 
simulated,  respectively,  the  buildup  of  water  in  the  lock  up  to  el  1+2.0 
and  development  of  uplift  pressures.  Uplift  pressures  equal  to  the 
hydrostatic  pressure  due  to  the  waterhead  in  the  lock  were  applied  at 
nodes  at  the  bases  of  monoliths  10-L  and  11-R  and  at  the  membrane  (see 
fig.  i+).  For  each  cycle,  three  iterations  were  performed  and  found  to 
be  satisfactory  for  a convergent  solution. 

Comparisons 

Settlements 

15.  Fig.  8 shows  a comparison  between  the  predicted  vertical 
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NOTE:  NUMBERS  INDICATE  FE  SEQUENCES  (TABLE  2) 


Fig.  8.  Comparisons  between  predicted  and  observed 
settlements , monolith  10-L 

displacements  measured  at  nodes  near  typical  settlement  plates  in  the 
backfill  adjacent  to  monolith  10-L  (fig.  3)  and  the  observed  settlements 
of  the  plates.  The  FE  results  are  shown  for  four  stages:  completion 
of  lock  wall  construction,  completion  of  backfill,  filling  the  lock,  and 
after  the  application  of  uplift  pressures.  The  correlation  between 
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observations  and  predictions  is  considered  to  be  satisfactory.  The  ob- 
servations showed  that  after  water  was  introduced  into  the  lock,  the 
settlements  (up  to  1972)  remained  essentially  the  same  or  decreased  in 
magnitude.  The  decrease  may  be  due  to  the  effect  of  uplift  pressures. 
The  FE  computations  with  full  uplift  pressure  also  showed  a decrease  in 
predicted  settlements. 

16.  In  fig.  9a  are  shown  observed  and  computed  displacements  at 
various  locations  in  the  backfill  side  that  took  place  between  the 
period  lit  August  1968  (end  of  backfilling)  to  7 January  1972  (during 
operations).  The  comparisons  for  plate  2 (fig.  3)  below  the  heel  of 
the  monolith  are  not  shown  for  clarity;  the  computed  settlement  of 
this  plate  was  about  -0.0 it  ft  whereas  the  observed  value  was  on  the 
order  of  -0.055  ft.  Although  the  FE  computations  do  not  provide  exact 
simulation  of  natural  and  man-made  conditions , the  computed  vertical 
displacements  at  typical  nodes  in  the  backfill  for  a short  portion, 
from  completion  of  the  backfill  (August  1968)  to  filling  the  lock 
(November  1968),  of  the  foregoing  period  seem  to  follow  the  trends  of 
the  observations  and  are  of  the  same  order  of  magnitude,  fig.  9a.  These 
nodes  are  in  the  neighborhood  of  the  locations  of  the  corresponding 
settlement  plates  in  the  backfill,  fig.  3. 

17.  Fig.  9’o  shows  computed  vertical  displacement  patterns  for 
various  stages  of  construction  at  two  sections:  at  the  base  of  the 
structure  and  at  section  B-B  in  the  foundation.  Dewatering  (stage  l) 
showed  settlements  whereas  at  the  end  of  excavation  (stage  U),  the 
highest  rebound  occurred.  Subsequently  the  foundation  showed  settle- 
ments as  the  construction  proceeded.  After  the  structure  was  com- 
pleted, the  distribution  of  displacements  was  essentially  symmetrical. 

At  the  end  of  backfill  (stage  10 ),  the  foundation  zone  below  monolith 
10-L  settled  more  than  the  zone  below  monolith  11-R  because  the  backfill 
was  provided  only  behind  monolith  10-L.  Under  the  uplift  pressures 
(stage  12),  the  foundation  settlements  showed  a decrease.  The  final 
settlements  are,  however,  less  than  the  total  rebound  after  excavation. 
This  seems  reasonable  since  the  approximate  weight  of  material  exca- 
vated was  about  57»000  tons  whereas  the  total  applied  load,  including 
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the  structure,  backfill,  and  water,  was  about  35,000  tons.  The  weight 
of  the  excavated  material  and  backfill  in  these  computations  includes 
the  entire  area  that  was  affected  by  excavation  and  backfill,  as  shown 
in  the  FE  mesh,  fig.  4. 

18.  Computed  vertical  settlements  at  nodes  199  and  483,  fig.  9, 
beneath  monoliths  11-R  and  10-L,  respectively,  are  shown  in  fig.  10. 
The  average  observed  vertical  displacements  of  settlement  plates  1 and 
2,  monolith  10-L,  fig.  3,  after  completion  of  lock  construction  and 
backfill  shown  in  fig.  10  compare  favorably  with  the  corresponding  FE 
results,  sequences  7 and  10  in  table  2.  The  displacements  under  both 
monoliths  are  about  the  same  up  to  the  end  of  lock  buildup  (stage  7), 
with  monolith  10-L  showing  somewhat  larger  displacements.  This  may  be 


Fig.  10.  Settlement  versus  construction  sequences  at  typical 
nodes,  199  and  483,  fig.  9 
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due  to  the  fact  that  the  soil  beneath  10-L  contains  a greater  depth  of 
relatively  weaker  Tertiary  clay  as  compared  with  the  depth  of  substratum 
sands,  fig.  L,  than  does  the  soil  beneath  11-R.  Monolith  10-L  shows 
significantly  larger  displacements  after  the  backfill  behind  10-L  is 
completed. 

Loads  in  piles 

19.  Typical  results  for  piles  battered  out  (B  piles)  and  battered 

in  (A  piles)  are  shown  in  fig.  11.  Distributions  of  loads  along  each 

3 

of  the  instrumented  rows  as  computed  by  Hrennikoff ' s method  are  also 
shown  in  fig.  11. 

20.  An  approximate  value  of  load  in  a pile  was  obtained  on  the 
basis  of  the  computed  vertical  stress  in  the  element  at  the  butt  of 

the  pile.  This  stress  was  multiplied  by  the  width  (l  ft)  of  the  equiva- 
lent pile  and  length  (1*0  ft)  of  the  monolith  (fig.  6),  and  the  product 
was  divided  by  the  number  of  piles  in  the  row  along  the  length. 

21.  As  shown  in  fig.  11,  the  measurements  were  recorded  at  vari- 
ous times  after  the  structure  was  completed.  In  some  cases,  measure- 
ments for  only  one  pile  location  were  available.  The  computed  distribu- 
tions of  loads  are  in  good  agreement  with  the  measured  distributions  for 
the  cases  where  sufficient  data  were  available.  The  measurements  showed 
that  the  outer  row  carried  higher  loads  than  the  piles  in  the  center. 
This  observation  is  consistent  with  previous  reports  for  friction  pile 
groups.^*'’*'^  The  predictions  for  the  completion  of  backfill,  particu- 
larly for  monolith  10-L,  fig.  lib,  also  support  this  observation. 

22.  The  distributions  of  loads  given  by  Hrennikoff' s method  are 
always  linear  and  do  not  in  general  agree  with  the  observations.  This 
may  be  due  to  the  fact  that  Hrennikoff' s method  does  not  account  for 
nonlinear  behavior  of  soils  and  interfaces  and  the  interaction  behavior 
of  the  system. 

23.  Fig.  12a  shows  computed  distributions  of  stresses  in  the 
piles  under  monoliths  10-L  and  11-R  for  the  case  of  piles  battered  out. 
As  expected,  the  piles  on  the  backfill  side  show  higher  stresses  than 
the  piles  below  monolith  11-R.  The  stresses  increase  in  magnitude  up 
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load,  kips  load  kips 


to  sequence  11,  water  in  the  lock;  thereafter,  under  uplift  pressure 
the  stresses  in  the  piles  decrease. 

Soil  pressures 


i 

\ 


2h.  Computed  distributions  of  vertical  soil  stresses  in  the  ele- 
ments at  section  B-B  in  the  foundation  for  typical  sequences  are  shown 
in  fig.  12b.  The  soil  stresses  are  higher  under  monolith  10-L,  and 
they  increase  in  magnitude  up  to  sequence  11.  Under  uplift  pressure, 
the  soil' stresses  decrease  in  magnitude. 

25.  Computed  distributions  of  vertical  soil  stresses  in  the  ele- 
ments immediately  below  the  slab  of  the  lock  are  shown  in  fig.  12c. 
Lateral  pressure  in  backfill 

26.  Computed  distributions  of  the  vertical  and  horizontal 
stresses  in  the  elements  along  the  section  Y-Y  in  the  backfill  region 
are  shown  in  fig.  13.  The  vertical  stresses  follow  essentially  the 
linear  variation,  = yD  , for  a major  portion  of  the  depth,  where 

a = vertical  stress  , y = unit  weight  of  backfill  , and  D = depth  of 
overburden  . In  the  lower  portion,  however,  the  pressure  shows  a 
significant  decrease  as  compared  with  the  linear  distribution.  One  of 
the  reasons  for  such  a reduction  can  be  the  existence  of  arching  near 
the  base  of  the  wall.  The  magnitudes  of  decrease  by  small  amounts 

from  sequence  10  to  12  in  the  upper  zones,  but  increase  by  small  amounts 
in  the  lower  zone. 

27.  The  horizontal  stress,  , also  decreases  by  somewhat 

larger  amounts  from  sequence  10  to  12  in  the  upper  zones  at  constant 

elevation,  and  increases  by  small  amounts  in  the  lower  zone.  For  the 

c 

major  portion  of  the  upper  zone,  the  linear  distribution,  0 = K yD  , 

c 

where  K = coefficient  of  lateral  earth  pressure  for  the  backfill 
o 

(0.1+5)  , lies  in  the  vicinity  of  the  computed  distributions. 

Drag  forces 

28.  Differential  settlement  between  the  backfill  and  the  lock 

I* 

can  cause  drag  forces  which  can  increase  loading  on  the  foundation. 

The  drag  forces  were  computed  from  the  FE  results  by  summing  the  shear 
forces  induced  in  the  elements  along  the  section  Y-Y,  fig.  13,  adjacent 
to  the  heel  of  the  monolith  10-L.  The  net  value  of  such  drag  forces 


Distribution  of  computed  pressures  in  backfill 


after  the  end  of  the  sequence  11  (table  2)  was  obtained  as  equal  to  about 
10  tons  per  foot.  The  properties  of  the  backfill  were  chosen  as  0 
= ho  deg  and  = 0.1*5  (table  l). 

29.  The  net  load  applied  due  to  the  weight  of  monolith  10-L  and 

the  weight  of  the  backfill  over  the  heel  was  estimated  to  be  equal  to 

1* 

about  130  tons  per  foot.  Hence,  the  drag  force  as  a percentage  of 
the  total  load  appears  to  be  about  8 percent;  this  value  compares 
favorably  with  that  of  10  percent  reported  in  reference  h. 

Comments  on  Formulation  and  Codes 

g 

30.  The  computer  codes  developed  previously  and  used  for  this 
study  can  yield  satisfactory  results  for  stresses  and  deformations 

in  the  foundation  and  surrounding  soils.  The  formulation  and  codes, 
however,  may  involve  certain  difficulties.  The  stresses  and  deforma- 
tions in  the  structure  itself  and  at  some  locations  in  the  interface 
may  not  be  realistic;  consequently,  interaction  effects  may  not  be 
studied  precisely.  In  the  formulation,  it  is  required  to  employ  "air" 
elements  during  various  sequences  of  operations.  Excessive  displace- 
ments can  occur  at  the  nodes  of  these  elements  in  comparison  with  the 
elements  in  the  intact  (foundation)  elements.  This  may  be  due  to  the 
fact  that  the  relative  magnitudes  of  the  component  of  the  overall  stiff- 
ness matrix  pertaining  to  the  air  elements  are  much  smaller  in  compari- 
son with  the  remaining  portion  of  the  stiffness  matrix  related  to  the 
intact  elements.  As  a consequence,  some  ill-conditioning  may  enter  the 
system.  The  stresses  in  elements  situated  symmetrically  may  not  be 
equal,  and  the  constitutive  models  for  soils  and  interfaces  may  require 
improvements.  Some  of  these  aspects  are  being  examined  currently  at  WES. 
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PART  V:  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


31.  The  three-dimensional  lock  pile  foundation  system  for  the 
gravity-type  Columbia  Lock  was  idealized  as  a two-dimensional  problem. 
The  system  was  simulated  by  an  equivalent  plane-strain  idealization  in 
which  a row  of  H-piles  in  the  group  was  represented  by  a strip  of  unit 
width  and  length  with  equivalent  (axial)  stiffness  properties.  The  se- 
quences of  construction — in  situ  stresses,  dewatering,  excavation,  con- 
struction of  piles  and  lock,  backfilling,  filling  the  lock  with  water, 
and  uplift  pressure — were  simulated  in  the  incremental  FE  analysis. 
Nonlinear  deformation  characteristics  of  soils  and  interfaces  were  de- 
rived from  laboratory  triaxial  and  direct  shear  tests,  respectively, 
and  were  introduced  into  the  FE  procedure. 

32.  Even  with  the  foregoing  idealization  to  a two-dimensional 
system  and  approximations  to  the  soil  and  interface  parameters,  the 
predicted  behavior  in  terms  of  progressive  settlements  and  distribu- 
tions of  loads  in  the  piles  showed  satisfactory  correlation  with 
corresponding  observations.  The  distribution  of  loads  in  the  piles 
from  the  conventional  Hrennikoff's  method  did  not  show  good  agreement 
with  the  pattern  of  load  distributions,  whereas  the  FIE  method  showed, 
in  general,  improved  correlations.  It  should  be  noted  that  the  loading 
conditions  and  material  properties  assumed  in  the  design  with  Hrenni- 
koff's method  may  be  different  from  those  simulated  in  the  FE  analysis. 
The  comparison  herein  is  made  only  as  a qualitative  trend. 

33.  The  concept  of  substituting  for  the  three-dimensional  system 
an  equivalent  two-dimensional  system  for  the  FE  analysis  seems  to  be 
satisfactory  from  a practical  viewpoint.  It  shows  the  engineer  that 
satisfactory  predictions  can  be  obtained  by  using  this  idealization 
which  can  provide  significant  savings  in  terms  of  human  and  computa- 
tional efforts.  The  FE  method  can  provide  useful  insight  into  and  in- 
creased understanding  of  the  behavior  of  battered  pile  groups  for 
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foundations  of  gravity-type  structures , and  can  be  used  as  a tool  in 
analysis  and  design. 


Recommendat ions 


31*.  The  investigation  reported  herein  was  directed  essentially 
to  the  study  of  stresses  and  deformations  in  the  foundation  soils.  As 
remarked  earlier,  the  present  formulation  and  code  may  show  local  errors, 
and  the  preciseness  of  interface  stresses  may  not  be  satisfactory.  It 
is  proposed  that  the  formulation  and  code  be  examined  thoroughly  and  that 
possible  improvements  be  introduced.  The  mathematical  properties  of 
the  interface  element  may  require  close  scrutiny. 

35.  The  interaction  effects  between  the  structure  and  soils  are 
highly  relevant  from  viewpoints  of  the  design  engineer.  With  the  fore- 
going suggestion  in  paragraph  3^,  it  will  be  very  useful  to  adopt 
finer  meshes  for  the  structure  and  compute  detailed  distribution  of 
stresses  in  the  concrete  structure.  This  will  permit  computation  of 
bending  moments  and  shear  forces  at  critical  locations.  A comparison 
of  such  design  quantities  with  the  corresponding  quantities  based  on 
conventional  methods  can  be  relevant. 

36.  The  topic  of  group  effects  in  pile  groups  is  of  importance 
from  the  viewpoint  of  design  analysis.  It  is  proposed  that  the  finite 
element  procedure  be  extended  and  used  to  examine  the  group  effects 
and  distribution  of  loads  in  the  piles  in  a group  and  the  results  com- 
pared with  the  loading  capacity  of  such  single  piles. 

37.  A parametric  study  to  examine  effects  of  such  factors  as 

0 

different  values  of  Kq  , drained  and  undrained  parameters,  various 
mesh  layouts , and  interfaces  at  the  base  should  be  performed  to  assist 
the  user  in  the  proper  choice  of  these  quantities.  Moreover,  it  is 
possible  to  perform  trial-and-error  analyses  in  which  both  the 
Hrennikoff's  and  FE  methods  are  used  under  various  combinations  of 
such  factors  as  material  properties,  loadings,  and  geometry.  These 
can  lead  to  better  criteria  for  design  analysis  of  gravity  locks. 


LITERATURE  CITED 


1.  U.  S.  Array  Engineer  District,  Vicksburg,  "Columbia  Lock,  Pile 
Tests — Lock  Site,"  Design  Memorandum  No.  6,  Supplement  No.  2,  Dec 
1963,  Vicksburg,  Miss. 

2.  , "Columbia  Lock,  Masonry  and  Embedded  Metals,"  Design 

Memorandum  No.  6,  Mar  1963,  Vicksburg,  Miss. 

3.  Hrennikoff,  A.,  "Analysis  of  Pile  Foundations  with  Batter  Piles," 
Transactions,  American  Society  of  Civil  Engineers,  Vol  115,  Paper 
No.  2 1*01 , 1950,  pp  351-381. 

U.  Montgomery,  R.  L.  and  Sullivan,  A.  L.,  Jr.,  "Interim  Analysis  of 
Data  from  Instrumentation  Program,  Columbia  Lock,"  Miscellaneous 
Paper  S-72-30,  Aug  1972,  U.  S.  Army  Engineer  Waterways  Experiment 
Station,  CE,  Vicksburg,  Miss. 

5.  Sherman,  W.  C.,  Jr.,  "The  Behavior  of  Lock  Walls  Supported  on 
Batter  Piles,"  Proceedings,  8th  International  Conference  on  Soil 
Mechanics  and  Foundation  Engineering,  Moscow,  Oct  1973. 

6.  Worth,  N.  L.  et  al.,  "Pile  Tests,  Columbia  Lock  and  Dam,  Ouachita 
and  Black  Rivers,  Arkansas  and  Louisiana,"  Technical  Report  3-7^1, 
Sep  1966,  U.  S.  Army  Engineer  Waterways  Experiment  Station,  CE, 
Vicksburg,  Miss. 

7.  U.  S.  Army  Engineer  District,  Vicksburg,  "Ouachita-Black  Rivers 
Navigation  Project,  Ouachita  River,  Louisiana,  Columbia  Lock  and 
Dam,"  Periodic  Inspection  Report  No.  1,  Oct  1972,  Vicksburg,  Miss. 

8.  Clough,  G.  W.  and  Duncan,  J.  M. , "Finite  Element  Analyses  of  Port 
Allen  and  Old  River  Locks,"  Contract  Report  S-69-6,  Sep  1969,  U.  S. 
Army  Engineer  Waterways  Experiment  Station,  CE,  Vicksburg,  Miss.; 
prepared  under  Contract  No.  DACW39-68-C-00U0  by  Lower  Mississippi 
Valley  Division,  Vicksburg,  Miss. 

9.  Desai,  C.  S.,  "Finite  Element  Method  for  Design  Analysis  of  Deep 
Pile  Foundations"  (in  preparation),  U.  S.  Army  Engineer  Waterways 
Experiment  Station,  CE,  Vicksburg,  Miss. 

10.  , "Nonlinear  Analyses  Using  Spline  Functions,"  Journal, 

Soil  Mechanics  and  Foundations  Division,  American  Society  of  Civil 
Engineers,  Vol  97»  No.  SM10,  Oct  1971,  PP  1^61-1^80. 

11.  , "Numerical  Design  Analysis  for  Piles  in  Sands,"  sub- 

mitted to  be  published  in  Journal,  Soil  Mechanics  and  Foundations 
Division,  American  Society  of  Civil  Engineers. 

12.  Desai,  C.  S.  and  Abel,  J.  F. , Introduction  to  the  Finite  Element 
Method,  Van  Nostrand  Reinhold  Co.,  New  York,  1972. 

13.  Desai,  C.  S.  and  Holloway,  D.  M. , "Load-Deformation  Analysis  of 
Deep  Pile  Foundations,"  Proceedings,  Symposium  on  Applications  of 
the  Finite  Element  Method  in  Geotechnical  Engineering,  C.  S.  Desai, 
ed. , Sep  1972,  pp  629-656. 


28 


Duncan,  J.  M.  and  Chang,  C.  Y.,  "Nonlinear  Analysis  of  Stress  and 
Strain  in  Soils,"  Journal,  Soil  Mechanics  and  Foundations  Division, 
American  Society  of  Civil  Engineers,  Vol  9^,  No.  SM5,  Sep  1970, 
pp  1629-1^53. 

Kulhawy,  F.  H. , Duncan,  J.  M. , and  Seed,  H.  B.,  "Finite  Element 
Analyses  of  Stresses  in  Movements  in  Embankments  During  Construc- 
tion," Contract  Report  S-69-8,  Nov  1969,  U.  S.  Army  Engineer  Water- 
ways Experiment  Station,  CE,  Vicksburg,  Miss.;  prepared  under 
Contract  No.  DACW39-68-C-0078,  by  Soil  Mechanics  and  Bituminous 
Materials  Laboratory,  University  of  California,  Office  of  Research 
Services,  Berkeley,  Calif.,  Report  No.  TE-69-k. 

Chellis,  R.  D. , Pile  Foundations,  2d  ed. , McGraw-Hill,  New  York, 

1961. 


Table  1 


Hyperbolic  Parameters  and  Properties  for  Different  Materials 


Material* 

V 

Vf 

y , lb/ft3 

Rf 

0 0 

0 

C 

K 

K 

u 

n 

1 

0.10 

0.10 

— 

— 

— 

— 

— 

— 

— 

0.0 

2 

0.30 

0.1+8 

120.0 

0.90 

0.70 

26 

0 

1000 

2000 

0.0 

3 

0.30 

o.i+8 

58.0 

0.90 

0.70 

26 

0 

1000 

2000 

0.0 

1* 

0.30 

0.1+8 

53.0 

0.85 

0.1+5 

l+o 

0 

1160 

1750 

0.5 

5 

0.30 

0.1+8 

58.0 

0.85 

1.10 

30 

0 

1+00 

800 

0.5 

6 

0.28 

0.28 

61.0 

1.00 

1.00 

— 

— 

— 

— 

— 

7 

0.28 

0.28 

66.0 

1.00 

1.00 

— 

— 

— 

— 

— 

8 

0.10 

0.10 

0.0 

— 

— 

— 

— 

— 

— 

0.0 

9 

0.20 

0.20 

150.0 

— 

1.00 

— 

— 

— 

— 

— 

10 

0.30 

0.1+8 

115.0 

0.85 

0.1+5 

1+0 

0 

580 

860 

C.5 

11 

0.U8 

0.1+8 

62.1+ 

1.00 

1.00 









0.0 

Note:  v = Poisson's  ratio. 

= Poisson's  ratio  at  failure. 
y = unit  weight. 

R^,  = failure  ratio. 

= coefficient  of  earth  pressure  in  compression. 

0  = Mohr-Coulomb  (angle  of  internal  friction)  parameter, 
c = Mohr-Coulomb  cohesive  strength  parameter. 

K = hyperbolic  loading  parameter. 

Ku  = hyperbolic  unloading  parameter, 
n = experimentally  determined  parameter. 

* Material  Identification 


1 Air 

2 Backswamp  clay  above  water  table 

3 Backswamp  clay  below  water  table 

U Substratum  sand 

5 Tertiary  clay 

6 Equivalent  piles  - 11-R 

7 Equivalent  piles  - 10-L 

8 Air  as  replacement 

9 Concrete 

10  Backfill  sand 

11  Water 


Table  2 

Sequences  Simulated  in  Finite  Element  Analyses 


Operation 

Details 

Sequences 

Initial  stresses 



__ 

Dewatering 

From  el  3^.0*  to  el  5.0 

1 

Excavation 

In  three  stages : 

El  58.0  to  1+2.0 

2 

El  1+2.0  to  26.0 

3 

El  26.0  to  8.0 

1+ 

Placement  of  piles 

— 

5a 

Lock  construction 

In  three  stages : 

El  8.0  to  26.0 

5b 

El  26.0  to  1+2.0 

6 

El  1+2.0  to  61+.0 

7 

Backfill 

In  three  stages : 

El  8.0  to  26.0 

8 

El  26.0  to  1+2.0 

9 

El  1+2.0  to  61+.0 

10 

Filling  of  the  lock  — 11 


Development  of  uplift  pressure 


12 
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In  accordance  with  ER  70-2-3,  paragraph  6c(l)(b), 
dated  15  February  1973,  a facsimile  catalog  card 
In  Library  of  Congress  format  is  reproduced  below 

K' 


Desai,  Chandrakant  S 

Finite  element  analysis  of  the  Columbia  Lock  pile 
foundation  system,  by  C.  S.  Desai,  L.  D.  Johnson  ..-and,' 

C.  M.  Hargett.  Vicksburg,  U.  S.  Army  Engineer  Waterways 
Experiment  Station,  1974. 

1 v.  (various  pagings)  illus.  27  cm.  (U.  S.  Water- 
ways Experiment  Station.  Technical  report  S-74-6) 

Sponsored  by  U.  S.  Army  Engineer  District,  Vicksburg. 
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